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We present a general theory of mixing for an arbitrary number of fields with integer or half-integer 
spin. The time dynamics of the interacting fields is solved and the Fock space for interacting fields 
is explicitly constructed. The unitary inequivalence of the Fock space of base (unmixed) eigenstates 
and the physical mixed eigenstates is shown by a straightforward algebraic method for any number of 
flavors in boson or fermion statistics. The oscillation formulas based on the nonperturbative vacuum 
are derived in a unified general formulation and then applied to both two and three flavor cases. 
Especially, the mixing of spin-1 (vector) mesons and the CKM mixing phenomena in the Standard 
Model are discussed emphasizing the nonperturbative vacuum effect in quantum field theory. 



I. INTRODUCTION 



O 
O 

<N 

, The mixing of quantum fields plays an important role in the phenomenology of high-energy physics 0,O,E8|. 
qq • Mixings of both K°K° and B°B° bosons provide the evidence of CP violation in the weak interaction jl4| and 7777/ 

(N '■ 

boson mixing in the SU(3) flavor group gives a unique opportunity to investigate the nontrivial QCD vacuum and 

(N 

fill the gap between QCD and the constituent quark model. In the fermion sector, neutrino mixing and oscillations 
are the likely resolution of the famous solar neutrino puzzle p6| , |r7| , p6| . In addition, the standard model incorporates 
the mixing of fermion fields through the Kobayashi-Maskawa (CKM) mixing of three quark flavors, a generaliztion of 
, the original Cabibbo mixing matrix between the d and s quarks [jll |22 24|,|2{|. Therefore, careful theoretical analyses 

o ■ 

of the mixing problem in quantum field theory is an important step toward understanding the many-body aspects of 
Oh, high-energy phenomena and their relationship to other areas of physics involving phase transitions. 

Moreover, the theory of mixing fields touches important, yet not fully answered, fundamental question about the 
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of the interacting fields and the free fields. This is different from the perturbation theories where the vacuum state of 



quantization of the interacting fields. The mixing transformation introduces very non-trivial relationships between the 
interacting and non-interacting (free) fields, which lead to a unitary inequivalence between the two Fock spaces 



interacting fields is equal to the vacuum of free fields up to a less essential phase factor e lSo |l0|]r2] , |2l| . The mixing 
of quantum fields is one of the cases that can be solved nonperturbatively in the quantum field theory. Thus, it also 
allows to investigate the accuracy of perturbation theory. For instance, the dynamics of a mixed-field Hamiltonian 
can be used for a partial summation of regular perturbation series as well as an improvement of the accuracy in 
perturbation theory. 

Recently, importance of the mixing transformations has prompted a fundamental examination of them from a 
quantum field theoretic perspective. The investigation of two-field unitary mixing in the fermion case demonstrated 
a rich structure of the interacting-field vacuum as SU(2) coherent state and altered the oscillation formula including 
the antiparticle degrees of freedom. Momentum dependence of mixing, existence of correlated antiparticle beam 
and additional high-frequency oscillation terms have been found and at the same time the vacuum condensates have 



1 



been analyzed for fermions P,p|,p|.p| Jl9| . |20| | . Subsequent analyses for the boson case revealed similar features but much 
more complicated vacuum structure for interacting fields [p|j7|p^|. Especially, the pole structure in the inner product 
between the mass vacuum and the flavor vacuum was found and related to the convergence limit of perturbation series 
p3[ . Attempts to look at the mixings of three- fermion case have also been carried out Q,^,^] . 

In this paper, we extend the previous analyses of mixing phenomena and work out a unified theoretical framework 
for an arbitrary number of flavors with any integer (bosons) or half-integer (fermions) spin statistics. We build the 
representation of mixing transformation in the Fock space of quantum fields and demonstrate how this can be used 
to obtain exact oscillation effects. We then use the developed framework to carry out calculations of two-field and 
three-field unitary mixings for the typical spin (i.e. 0,1/2 and 1) cases. We also comment on the use of mixed-field 
solution to improve the perturbation series of mixing effects. 

The paper is organized as follows. In Section II, we define the ladder operators for flavor fields and explicitly show 
the unitary inequivalence between the flavor Fock space and the Fock space of mass-eigenstates. In Section III, we 
find the time dynamics of the flavor ladder operators and derive general expressions for the particle condensations and 
the number operators as functions of time. We also present some remarks on Green function method in the mixing 
problem. We then specifically consider, in Section IV, the mixing of two spin-1 fields (vector mesons) along with the 
mixing of spin-1/2 fields and show the consistency with previously known results. Summary and conclusion follow in 
sections V. In Appendix A, the mixing parameters are shown explicitly for the spin 0, 1/2 and 1. In Appendix B, we 
present a derivation of the flavor vacuum state by solving an infinite system of coupled equations which appears as a 
condition of the vacuum annihilation. In Appendix C, we summarize our results of the three-field mixing for the spin 
0,1/2 and 1 using the SU(3) Wolfcnstcin parametrization. 

II. THE THEORY OF QUANTUM FIELD MIXINGS 

In this section, we consider the mixing problem for N fields of fermions or bosons. To discuss the dynamics of the 
flavor (mixed) fields, we define a flavor field </> M as a mixture of the free fields (fj (j = 1, 2, . . . , N); i.e. 
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where U^j is a unitary mixing matrix element. We use the latin indices i,j,k,... to label the fields of mass- 
eigenstates and the greek indices /i, v, £, . . . to label the flavor fields. We also denote and <p as the entire columns 
<f> = (0i,02 ■ ■ ■ ,0jv) T and (p — (<fi,(p2 ■ ■ ■ ,Pn) T , respectively. The evolution of the fields M is generated by the 
Hamiltonian of the form R 



1 When there is an additional interaction Hamiltonian for given by Hi = cfrWcf), the Hamiltonian of the system is, of course, 
extended to H{<j>) = H{<j>) + Hi = H (<j>) + <j> ] M<f> + ftW<j>. Then Hf ree (0) is also extended to Hf Tee {ip) + 0^U^WU(p. 
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H($) = H free {<p) = H free (U^) = H (4>) + ftM4>, (2.2) 

where Hf ree ((p) is the free field Hamiltonian for ipi with the corresponding mass eigenvalues m,, Ho(<f)) is the free 
flavor field Hamiltonian and M is a mixing matrix. 

The existence of the explicit relationship between free (ip) and flavor ((/>) fields, given by Eq.( |2.lD , allows us to work 
out the quantum-field theoretical solution to the problem given by 

j t ^=i[H$), (2.3) 

In fact, the solution of Eq.(fO) is contained in Eq.(|2.iP with the free field (ipi) given by 



where a.? (t) = e «-r and fc.g (t) = e 8e *£ 6.g with the standard equal time commutation/anticommutation 
relationships for bosons/fermions, i.e. 

K (t) , O*, (t)]± = [6a (t) , (t)]± = «a,a'. 



In Eq.(2.4), u\ and are the free particle and antiparticle amplitudes, respectively, and a is the helicity quantum 



number given by 

(n-s)ul a = <Tu}; a ,(n-s)vl a =avl a , (2.5) 

where s is the spin operator and n = k/\k\. We also define the following parameters that are useful in extracting the 
ladder operators from the field operators 



ka ko' — k — a —k—cr' 1 

ht i 8„„> = «SV - 



(2.6) 



ka a ' a ka —k — a' 



For the analysis of arbitrary flavor mass parametrizations, we use the general notation given by Eq.( |2.6| ) including 
both flavor and mass degrees of freedom. Although both indices \x and j are numbers running from 1 to N, the mass 
for the first index should be used as the flavor mass while the second index is for the mass eigenvalue m,j . One should 
note that H and h are both symmetric for bosons while H is symmetric and h is antisymmetric for fermions. The 
explicit representations of H and h are presented in the Appendix A for the spin 0,1/2 and 1 cases. 

Now, if A(Z7, t) is the representation of the mixing transformation defined in the equal time quantization, then 

<j>(t) = UCp{t) = A(U, tfy(t)A(U, t). (2.7) 

In the associate Fock-space, this corresponds to 

\a,t > f =A{U,t)' f \a,t > m (2.8) 
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where subscript / (to) is used to denote the flavor (mass) Fock-space. For the given time t, Eq.(E 1 2) can then be 
written as 



H{<t>{t)) = K{U,tyH free {(p{t))K{U,t). 



(2.9) 



As noticed from the two-field mixing analysis [pPj7|,^|, H{4>{t)) and H(<p(t)) cannot be in general related by the 
same operator at all times so that A(U, t) is essentially time dependent. The vacuum state of the flavor-fields, defined 
as the state with the minimum energy, is A(U, ty\0 > m and changes with time satisfying 



< a\H${t))\a > f = m < a\H(<p(t))\a > m > m < 0\H(<p(t))\0 > m =f< 0\H($(t))\0 >/ 



(2.10) 



We now define the ladder operators for the flavor fields as a =ig(T (t) = A(C/,t) t a^ (j (i)A(C/,i). Using linearity of the 
mixing transformation, we then can solve the explicit structure of a % a (t) without finding A(U, t) itself. 

Such approach in fact has been known for some time for the fermion case |2(J, where it was noticed that fermion 
ladder operators for spin 1/2 can be extracted from quantum fields by means of 



(2.11) 



Since the Fourier component tp-i{t) = ^ i_ fu- a.r (t) + v l ~ b\ - (t)\ is obviously a linear combination of 
tfi (x, t), one can express ladder operators as linear combinations of the initial fields. Using the linearity of Eq.([2.7|), 
we get 



/2e.. E 



W*) = uV{A{U,t)^~ k {t)A{U,t))^ 



(2.12) 



j ka 



For the bosons, however, the ladder operators are not separated as in the fermion case, e.g. 

1 



u 1 - ifi r(t) 



(a- r (t) + hS b 

V lfecr\ / k<T 



(2.13) 



and in general h l i ^ 0. Eq.(2.13) implies that particles and antiparticles in boson case can not be distinguished 



unless time dynamics is considered. To deal with this problem we define ladder operators for bosons by 



i — k — o 



(2.14) 



With Eqs.( 2.11 ) and ( 2.14 ), we then derive for fermionsf^ 



Here, we abbreviate the notations and b\ ^ as a.j and bU, respectively. Similar abbreviation is used for a M and b- 
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-W f^MULu . a . + /UMlLu .&t 



(2.15) 



= E 



/; ■ 



(2.16) 



and for bosons: 



/2e u 



a u = 



. V (u^u j efl+€j a + u A1 V e "~^ ^ ^ ^ - 



I1U.4-. 21. 



E 



2 h^U M bl 3 , , 



(2.17) 



ft-*. = 4^ E 4 , e -^a\ + «"t 

/ — 4- / — 



kcr' e v J —k — a — h — <7 £ m J / 



u, 



= E V-^V-" {H »j ru * jb _ j + 



1 _ / li. 



(2.18) 



Denoting the spin of the mixed fields as S, we can unify the expressions for both fermion and boson in an identical 
form as 



by defining 



«p = E ( a M a j + 0M b -j) . 

5-M = E(o^6-i + (-l) M /^at), 



(2.19) 



(2.20) 



where 



Tfjpr fcrmions, 



We also note from unitarity that 




bosons. 
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TJlWJr fcrmions 



(2.21) 



bosons. 



I KjT + = . fcrmions; 

I \ a njf ~ \Pm\ 2 = \ U m\ 2 . bosons 

so that one can treat a^j and /3 W - as cosine and sine for fermions (cosh and sinh for bosons), respectively: 

J cos(# w ) fermions 
"» ~ Um 1 coshf0..,O bosons ' 



wja\Unj) fermions 

icosh(^ w ) bosons ' 
sin(# w ) fermions 
sinh(6» w ) bosons 



(2.22) 



(2.23) 
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From the fact that Eqs.(2.15)-(2.18) serve as the mixing group representation, one can conclude that 



J M "HX ~ 3 3 



(2.24) 



regardless of m^. Using the formulas, presented in the Appendix A, this can be explicitly verified for S — 0, 1/2, 1 by 
calculating, for example, -g^-- In every case can be reduced to -g^- = 7^/(to m ), e.g. for fermions 



d9fxjf 36^j dm u dm 



dsin(9 MJ ,) dsin(g„) 

/(m M ) - /K) = 



drrifj, drrifj, cos(6> MJ /) cos(6' w ) 

so that 9 M =6^- 6j, where cos(^) = 2^(V e M + "V + V £ m ~~ "V) and sin (^) = 27^(V £ a< + "V - V e M ~ "V)- 

The introduced ladder operators are consistent with the representation of the mixing transformation in the Fock 
space: 

\ afl + l,t> f = al(t)\a fl ,t> f = A(U,tfal(t)A(U,t)A(U,tf\ ai ,t >m= A^i^K + 1,* >m, (2.25) 
and the flavor vacuum state satisfies 

5 M (t)|0,t >/= A([/,i)t a ^)A([/,i)A(f/,i)t|0 > m = 0. (2.26) 



While Eq.fl2J^) may be viewed as the result of expanding flavor fields 4>n{x) in the basis parametrized by free-field 



mass rrii, it was noticed that one may as well expand flavor fields in the basis with the flavor mass parameters m 



which correspond to choosing ut ,v ? as free-field amplitudes with the flavor mass (m (1 ) in Eqs.(2.11) and (2.14) 



In other words, for any A(U,t), A'(U,t) — 7(i) _1 A([7, t)I(t), that can be obtained by means of a similarity 
transformation mixing o^g ff (t) and fe^ - (i) but leaving their combination in <fr(k) unchanged (i.e. <fi^(k,t) = 



I(t) 1 4> l ^ l (k,t)I(t)), is also a representation of the mixing group. The ladder operators, defined by Eqs.( 2.15 )-( [2.1(j| ), 



therefore depend on the choice of I(t) or, equivalently, the "bare" mass m M assigned to the flavor fields which is called 
as a mass parametrization. 

Although there are different opinions about whether or not the measurable results of the theory depend on the 
mass parameters j^,|[[l9]|2^] , we note that the mass parametrization problem indeed is not specific to the quantum 
mixing, but can be revealed in the free field case as well as in the perturbation theory. As discussed in |p3| , when 
dealing with the free field problem defined by the free Hamiltonian 



one may still consider the change of the mass parametrization m — > m M defined in |?| by 



3(*) A _ T -i (f] ( « ^ m ( e^-^pt jto+t&Xs \ ( o(0) \ . . 



G 



where ik = \/k 2 + m^ 2 and = \Jk 2 + to 2 . Indeed, as we observe in p3], the number operator for the free fields in 
this transformation is not conserved, e.g. for fermions 

< N >= \{&,&Ht)}\ 2 = \\Pk\ 2 e- 1 ^ + |A fe |V £ *f , (2.29) 

that may lead to obviously wrong conclusion that the number of particles in the free field case is not observable 
quantity. 

This can be also understood mathematically, once we note that the above transformation is equivalent to the 
splitting of the initial hamiltonian into 

H Q = % + H'j= Jd 4 P ({ (p^(P^) ~ ro„ W} + (tV - m 2 ) , (2.30) 

where additional self-interaction term, responsible for oscillation of < N >, appears. Physically, the transformation, 



given by Eq.(2.28), can be viewed as a redefinition of the physical one-particle states. The tilde quantities correspond 
then to some new quasiparticle objects so that the tilde number operator describes a different type of particles and 
thus it doesn't have to be invariant under such transformation. Nevertheless, the charge quantum number is still 



conserved in the transformation, given by Eq.( 2.28| ). The situation here may be analogous to the representation of 



physical observables under the change of coordinate systems. Although the Casimir operator (e.g. S 2 in the spin 
observables) must be independent from the coordinate system, other physical operators (e.g. S x S y and S z ) do depend 
on the coordinate system. To compare the eigenvalue of S z between theory and experiment, one should first fix the 
coordinate system. Similarly, we think specific mass parameters should be selected from the physical reasoning to 
compare theoretical results (e.g. the occupation number expectation) with experiments. 

From the above example it is clear that the same mass parametrization problem is also present in the regular 



perturbation theory once one attempts to redefine the physical one-particle states as shown in Eq.( 2.28| ). Indeed, in 
the free threory and the perturbation theory this issue is resolved by the presence of the mass scale of well defined 
asymptotic physical states, which therefore fix the mass parameters. In the mixing problem, however, at least two 
feasible mass scales may be suggested either by the mass scale of the energy-eigenstates or by the flavor mass scale 
which corresponds to no self-interaction term in the hamiltonian, given by Eq.( |2.2[ ), and thus further discussion of 
this issue in the mixing problem is clearly necessary. We think the mass eigenvalues that can be measured from the 
experiments may be the natural choice for the mass scale in the given physical system. 

In any case, all the above unified formulation for any number of fields with integer or half-integer spin holds for the 



arbitrary mass parameter when ej = ^/k 2 + mf and e M = J k 2 + to 2 in Eqs.( 2.19 )-( 2.21 ) are understood as the 
energies of the free field tpi and the flavor field respectively. 

In the rest of this section, let us consider the explicit form of the flavor vacuum state. We obtain its structure by 
solving directly the infinite set of equations 

a„\0 > f = 0, K\0> f = 0. (2.31) 



We can express the flavor vacuum state as a linear combination of the mass eigenstates, i.e. in the most general form, 
l°>/= E n ,J ..„», *(«)(«) ■ ■ ■ (4)" (bU) h ■ ■ ■ (bl k ) h 10 > m , (2.32) 



(n),(l) 

with (n) — (niriina . . .). After applying Eq.(p3l|) to Eq.(|3|) we get an infinite set of equations given by 

^{a M B {nj+m) + PmB^.^x)) = 0, for all sets of (n), (I), (2.33) 
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where (rij ± 1) = (ni7i2 • • • fij ± 1 . . .). The solution of this problem is presented in the Appendix B. For the flavor 
vacuum state we find 

1 N 

|0 > f = - exp( ZiAb-jW >m, (2.34) 

i,3=l 

where Zij is an element of the matrix Z = —dr x fi. The normalization constant Z is fixed by / < 1 >/= 1; 
Z = dct 1/2 (l + ZzA for fermions and Z = det~ 1/2 (f - ZzA for bosons. The flavor Fock-space is then built by 
applying the flavor-field creation operators (aL6£) to the vacuum state |0 >/. 

We see that the flavor vacuum state has a rich coherent structure. This situation is different from the perturbative 
quantum field theory, where the adiabatic enabling of interaction is present and |0 > interacting™ |0 >/ ree . The 
nonperturbative vacuum solution renders non-trivial effects in the flavor dynamics as we will show in Section III. In 
particular, the normalization constant Z is always greater than 1 so that in the infinite volume limit, when density 
of states is going to infinity we have 

Z t ot = exp [ ^ f dkhx{Zz) J oo. (2.35) 



v (2?r) 3 , 

Thus, any possible state for the flavor vacuum shall have infinite norm in the free-field Fock space and therefore the 
flavor vacuum state cannot be found in original Fock space. The unitary inequivalence of the flavor Fock space and 
the original Fock space is therefore established, i.e. f < 0|0 > m ~ ^77 ~~ * The effect is essentially due to an infinite 
number of momentum degrees of freedom, which is analogous to the existence of phase transition in the infinite volume 
limit. 



III. TIME DYNAMICS OF THE MIXED QUANTUM FIELDS 

Now we have a closer look at the dynamics of quantum fields respresented by the ladder operators shown in 
Eq.( ^.lSj ). First of all, we note that only a ? and b i _£_ a operators and their conjugates are mixed together. We 
denote the set of quantum fields formed by all linear combinations of these operators and their products (algebra on 
a t , b i _j:_ a and h.c.) as a cluster fig with a particular momentum k and a particular helicity a. It follows that 
S7j ct 's are invariant under mixing transformation A(U, t) and we thus can treat each cluster independently from each 
other. 
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The time dynamics of the flavor fields is determined by the non-equal time commutation/anticommutation relation- 



ships for boson/fermion fields that can be derived from Eq.(2.19) using the standard commutation/anticommutation 
relationships for the original ladder operators; 



k,k' \ 

k 

G^(t) = [&_„(t),a„]± - £ (a; k f] uk , h-ke-^blJ + (-1) 2S 



(3.1) 



fc.fe' 



«k> 



The two matrices -F 1 and G represent the only nontrivial commutators/anticommutators in the sense that all others 
are either zero or can be written in terms of the elements of these matrices. It is useful to note that, for t = 0, Eq.([0]) 
shall be reduced to F MI/ (0) = and G M „(0) = 0. We also note that 



G^ u {t)* = —G v/ j,(t). 



(3.2) 



Eq. ( \A ) allows us to compute many mixing quantities directly. The time dynamics of the flavor-field ladder operators 
can be derived by writing them as a M (t) = E(//f ^"(0) + 9pvb— 1/ (0) + . . .). Then, one can get straightforwardly 
f* — [a u (0),al (t)]± — F v ^{—t) and — [b_ v (0),a,j, (t)]± = G vfl {—t) while all other coefficients are zeros: 

(*) = E (F^{t)a v + G ult (-t)bt v ) ; 



b-» (*) = E {F vlt {t)b-u + (-i) 2S G^(t)at) 



(3.3) 



We now consider the condensate densities of the definite-mass particles in the flavor vacuum [Z[ =/< 
0\a] (t)cii (t) |0 >/), the number of definite-flavor particles in the flavor vacuum {Z v = f < 0|4 (t) ~a v (t)|0 >/) and the 
particle number average for a single definite-flavor particle initial state, which is related in the Heisenberg picture to 
N pva = M < p\al (t) a v (t) \a > M , N pva = M < p\V_ v (t) 6-^ (i) |<r > jU . 

The free-field particle condensates in the flavor vacuum state are computed from the explicit form of the ladder 



operators given by Eq. ( 2.1£ ) as 



(3.4) 



In the following, the particle-antiparticle symmetry should be accounted for, so that a corresponding antiparticle 
quantity can be found from the particle expression after a necessary substitution (particles— »antiparticles and vice 
versa). Thus, the antiparticle condensate is given by the same quantity in Eq.(|3.4D. The definite-flavor particle 



condensates in the free-field vacuum is also given by Eq.(p.4|). 

Using Eq.(jDi), we get the flavor-field condensates in the flavor vacuum (Z v ) as 



Z v {t) =Y,\Gu^-tf 



(3.5) 



It is remarkable that this number is not zero but oscillates, demonstrating the oscillations of definite-flavor particles 

in the flavor vacuum. This effect reveals the unitary inequivalence of the flavor Fock-spaces for different times due to 

the time-dynamics of the flavor vacuum. 

The evolution of the particle (N pva ) and antiparticle (N pva ) number with flavor v can be found using the standard 

technique of normal ordering, i.e. moving annihilation operators to the right side and creation operators to the left 

side of the expression. With this technique, we obtain 

N pva (t) = [a p ,at (t)]±[a v (t) , a+]± + 5 pa < Q\a) v (t) a v (i) |0 >= 
= K P {t)F vtr (t)+6 pa Z v (t), 
N P ua(t) = (-1) 2S [a p , b-u (t)]±[bl v (t) , 5t]± + 8 pa < Ojfot (*) ~ hv (t) | >= ^ 
= (-1) 2S ' G vp {t)G va {ty + S pa Z u {t). 



The flavor charge Q pva = N pva — N pva [|7j— 19|| is then given by: 

Q pva = N pvc - N piJa = F* p {t)F va {t) - (-lf s G vp {t)G va (ty. (3.7) 

For a specific case of the number evolution in the beam with a fixed 3-momentum, we find: 

N pvp =< 0\a p at (t) a v (t) ajjo >= |F„ p (<)| 2 + Z u {t), 
N pvp =< 0\a p bl v (t) 6_„ (t) 5t|0 >= (-1) 25 \G„ p (t)\ 2 + Z u (t), (3.8) 
Q p „ P = \F up {t)\ 2 -(-l) 2S \G vp (t)\ 2 . 

We note that N pvp s as well as Q pvp s are in general dependent on the choice of mass parameter m M . 



We may explicitly see this in the example of the charge operator. According to Eq^^), we get 

Qpvp = EKfc<^ M - (-fMi e " i,t, )(^ Q * e " iet '' - (-^/^/We^)- 

-(-1) 2S E K.^e-^* - fta^e^'X^^e^' - /W^e""^) = 

k,k' 

(-1) 2S e^'+^ia^at^,^ - K k a pk a vk , p; k ,) (3 . 9) 
(-lf s e^'+^a^iat^, - f$ vk a vk ,) 

k,k' 

(-1) 2S (fi k <*uk> - atMie-^'+^Ppka;,, - e^'+^a^). 



Taking into account Eq.(2.23), we can write, e.g. for fermions (S=l/2) 

ot„k>o£ k + u k>0tk = U vk 'U* k (cos(6 vk >) cos(6 vk ) + sin(<W) sin(6U)) 
= U vkl U* k cos{6 vk > - 6 vk ) = U uk >U* k cos(6 kkl ), 

PZ k otuk> ~ o? vk p vk > = U uk ,U* k {cos{6 uk ,) sin(6U) - cos(9 u k')sui(6 vk )) 
= U v yU* k sm(6 vk - vk >) = U uk ,U* k sin(0 k > k ). 



Thus, wc find 



Qpvp = E U vk >Ul k U pk U* k , (cos 2 (6 kk > ) cos(uj k , k t) +icos(0fe' fe )cos(6' Alfe + pk ,) sm(u kk rt)+ 

k,k> (3.10) 
sm 2 (8 k , k )cos(n k , k t) - ism(9 k , k ) sm(9 pk + 9 pk >) sm(n kk ,t)); 
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where Q^j = ej + tj and = 6j — ej . This can be rewritten as 
Q^m = E Re(U„ k ,UZ k U^ k U* k ,)(cos 2 (6 kk ,)cos(u k > k t) - (-l) 2S sm 2 (8 k , k ) cos(n k , k t))+ 

k,k' /g -qS 

+ J] Im(Uv k ,U* k U^ k U* k ,){cos(9 kk ,) cos(9^ lk + 9^ k ,) sin(uj k , k t) - (-1) 2S sin(0 fe ' fe ) sin(f (tlfe + 0^) sin(fifc/ k t)). 
fe.fe' 



This formula is also valid for bosons with the substitution of cos — ► cosh, sin — * sinh. 

We see now that Q^ VI1 does not depend on the mass parameters only for real mixing matrices U^k f§|,[l9| . Otherwise, 
there is a nontrivial mass dependence from the imaginary part of U . Interestingly, even in the latter case, there is no 
dependence on the mass of the flavor field v (m„) but only on the mass of the initial flavor state /i. 



We also note that Eq.(B.S) may be viewed as a superposition of the two terms: p — > v propagation and back- 
ground vacuum contribution Z v . Thus, one may introduce the particle-particle and particle-antiparticle propagation 
amplitudes, respectively, defined by 

7V^(M) = [a„(t),at(o)]± = F vp {t) 

V p ^-*(k,t) = [b_ v (t),a p (0)] ± =G„ p (t). ^ 1 

Indeed, such propagation amplitudes appear from the flavor-field Green function f < 0(t = 0)\<f> v (k,t)<f>l(k,0)\0(t = 
0) > f for t > 0. Propagation functions, defined in this way, are clearly the Green functions of the mixed-field problem 
and obey the causality features relevant to such Green functions []. 

IV. TWO-FIELD UNITARY MIXING 
A. Vector Meson Mixing (S=l) 

We now consider the unitary mixing of 2 fields with spin 1 (vector mesons). U(2) parametrization consists of 4 
parameters: 3 phases that can be absorbed in the phase redefinition of fields and one essential real angle that is left, 
so that 

rr _ ( c °s(#) sin(0) \ . . 



-sin(0) cos(6>) 

Using Appendix A, we then define 7* = | ^y^± ^/f^) f° r c = ±1 and 

+ = 1 ^Ej-k 2 

- = 1 fc 2 +e^i 

for (7 = 0. For the free field mass rrii basis, 7^ = 7^ = 7+, 7^ = —721 = 7-- We use this basis in Section IV. 
The ladder mixing matrices a and (3 are given by 



(4.2) 



3 See Refs. HP-SI] for the discussion of the Green functions in the quantum theory of the mixed fields. 
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(4.3) 



cos (8) 7+ sin (8) 
— 7 + sin (9) cos (9) 

7_ sin (8) 
-7_ sin (0) 

For the flavor charge oscillation, we then obtain the result that is not dependent on the mass parametrization: 

Qui = 1 + sin 2 (28) ( 7 £ sin 2 - 7+ sin 2 (^)) , 

Q 121 = sin 2 (29) (<# sin 2 (^) - 7 2 sin 2 (^)) . ^ 

We see that this result, with an exception of greater complexity of j±, is identical to the case of spin J7|.p3| 



According to the above theory, in fact, this should be the case for the two-field mixing with any integer spin. For 
5 = 1 we see that an essential difference from the scalar/pseudoscalar meson mixing, such as the complication of 
momentum dependence of j± , occurs only for the mixing of longitudinally polarized particles. The mixing of transverse 
components is essentially same as in the case of spin-zero particles. 
The details of non-equal time commutators are given by 



F 



e -ieit cog 2 (0) + e -i £ 2t 7 2 sm 2 (Qj _ ei e 2 t 7 2 sin 2 ^ . J+ gin ^ CQS ^ (g-i«t „ e -ieit) 1 

7+ sin (9) cos (9) (e"" 2 * - e~ i£lt ) ; e~ ie2t cos 2 (9) + e"" 1 ^ 2 sin 2 (9) - e ieit -fl sin 2 (9) J ' ' 



7+7_ sin 2 (0) (e"" 2 * - e 162 ') 7_ sin (9) cos (6>) (e" lei * - - lt2t 



G ^ 7 _ sin (6») cos (6») (e" ie2t - e i£lt ) 7+7 _ sin 2 (f5) (e lei * - e~ leit ) J' ( ' 4,6 ^ 
The condensates of free-field particles are 

Z[ = Z' 2 = 7 2 sin 2 (0) (4.7) 

and the condensates of the flavor particles in the vacuum are 

Z x = 4 7 2 sin 2 (9) (cos 2 (9) sin 2 (2ui) + 7 2 sin 2 (0) sin 2 (S^)) , 

Z 2 = 4 7 2 sin 2 (9) (cos 2 (0) sin 2 (Oui) + 7 2 sin 2 (0) sin 2 (M)j . ^> 



The flavor vacuum structure is defined by the matrix Z: 



1 / -7+7- sin 2 (0) 7 _ cos (9) sin (0) 



(4.9) 



(cos 2 (6>) + 7 2 sin 2 (<?)) V 7- cos (8) sin (8) 7+7 _ sin 2 (0) 

with the normalization constant being Z = (1 — cos 2(g)+ 7 2 S i n 2 (g) ) -1 =1 + 7- sin 2 (0). 
The time evolution of the flavor particle number (if #1 was emitted) is given by: 

N ni = 1 + sin 2 (8) {8 7 2 cos 2 (0) sin 2 - 4 7 2 cos 2 (0) sin 2 (^) + 

+8 7 2 7 2 sin 2 (0)sin 2 (^)}, " (4.10) 

JV m = 4 7 2 sin 2 (8) (2 7 2 sin 2 (6) sin 2 (^f) + cos 2 (0) sin 2 (^)) , 

iViai = sin 2 (6») {4 7 2 cos 2 (0) sin 2 (^1) + 4 7 £ cos 2 (8) sin 2 (^) + 

+4 7 2 7 2 sin 2 (8) sin 2 (OuS) , " (4.11) 

JV121 = 4 7 2 sin 2 (9) (2 cos 2 (8) cos 2 (%^) + 7+ sin 2 (0) sin 2 (2$i)) . 

Also we note that the scalar and pseudoscalar case follows immediately from the above presentation when 7^ = 
i I . ± . J . In this respect, the spin-zero mixing is equivalent to the mixing of transverse components of vector 



fields, described by Eqs.(4.4), (4.5), (4.9) and ( 4.1C ). These results are consistent with the previously known results 
23l. 
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B. Fermion mixings (S=l/2) 



We also present here the calculations for s = 1/2 case. For the consistent notation with the previous works 0, 
we define 



jj _ \J (ei+m 1 )(e2+m 2 )+y/ (ei— mi)(e 2 -m 2 ) 

~~ 2 V £ i £ 2 ; (4.12) 

* = 17 7 /- . • 



The charge fluctuations are then given by 



Q m = 1 - sin 2 (20) (U 2 sin 2 (^) + V 2 sin 2 (^)) 



/m = sin 2 ' 



(20)(C/ 2 sin 2 (^)+y 2 sin 2 (V)) 
and the ladder mixing matrices are 



cos(0) C/sin(0) 
-E/sin(0) cos(0) 
/ Fsin(0)\ 
p sin (6») J 



(4.14) 

o J- 

which are same with the previously known results [jl|]6| . 

We can also give more details on the fermion mixing dynamics. The non-equal time anticommutators are given by 



g-ieit cog 2 (0) + e -ie 2 tjj2 gm 2 ^ + gie 2 ty2 ^2 ^ . g j n ^ CQg ^ ^-ie a t _ g-teit) 

f7 sin (0) cos (0) (e"" 2 * - e" 161 *) ; e" 162 * cos 2 (0) + e^ £l *[/ 2 sin 2 (0) + e" 1 *^ 2 sin 2 (0) 



f UV sin 2 (0) (e- ,e2 * - e" 2 *) V sin (0) cos (0) (e~ ieit - e" 2 *) \ , , 

\ V sin (0) cos (0) (e-" 2t - e ieit ) UV sin 2 (0) (e 4£lt - e" 161 *) J ' ^ ' 

The condensates of the free-field particles are 

Z[ = Z' 2 = V 2 sin 2 (0) (4.17) 

and the condensates of the flavor particles are 

Zx = W 2 sin 2 (0) (cos 2 (0) sin 2 (Oa*) + [/2 sin 2 (#) sin 2 (DM)) ; 
Z 2 = 4^ 2 sin 2 (0) (cos 2 (0) sin 2 (Oui) + C/ 2 sin 2 (0) sin 2 (^)) . 

The vacuum structure is defined by the matrix Z: 



cos 



2( 



g = -L ( -UVsm 2 (9) l/cos(0)sin(0) 



(0) + C/ 2 sin 2 (0) V^cos(0)sin(0) UVsin 2 (e) 



with the normalization constant being Z = — <>,„, , L> . ■>,„, = t — ..a 1 . ,,,„. . 

The time evolution of the flavor particle number (if #1 was emitted) is then given by: 

Nni = 1-4U 2 sin 2 (0) cos 2 (0) sin 2 (^f) , 
N ln = W 2 sin 2 (0) cos 2 (0) sin 2 (%^) , 



(4.19) 



Nth = 4 sin 2 (0) {U 2 cos 2 (0) sin 2 + V 2 cos 2 (0) sin 2 (Ou*) + 

+C/ 2 F 2 sin 2 (0)sin 2 (%i)}, " (4.20) 

iV 121 = 4(7 2 ^ 2 sin 4 (0) sin 2 (Oul) . 



4 In our notation, [/ = 7+, V = 7_ 
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V. CONCLUSION 



The quantum field mixing effects may be understood by considering interplay between the two Fock-spaces of 
the free-fields and the interacting fields. As demonstrated in the 2-field mixing treatment, this interplay is highly 
non-trivial and gives rise to a deviation from the simple quantum mechanical approach due to the high-frequency 
oscillations and the antiparticlc component in the system. 

We have now extended the previous results and presented a solution without approximations for the quantum field 
theory of mixings in the arbitrary number of fields with boson or fermion statistics. As one might have expected from 
the previous 2-field treatment ||[||7],|23| , all results fall into the same scheme and can be easily unified. We investigated 
the field time dynamics by calculating unequal time commutators and discussed the propagation functions. We found 
an explicit solution for the interacting field Fock space and the corresponing vacuum structure that turned out to 
be a generalized coherent state. We then showed the unitary inequivalence between the mixed-field Fock space and 
the free-field Fock space in the infinite volume limit. After we built a formal calculational framework, we applied it 
to solve mixing dynamics of 2 vector mesons(S f = 1) and fermions(S' = 1/2). We found that the scalar/pseudoscalar 
(S = 0) boson mixing is the same as the mixing of transverse components of the vector fields, while for the longitudinal 
component of the vector field we found richer momentum dependence than in the spin-zero case. 

However, from the application of our approach to 3-fermion/boson mixing cases, which we summarize in Appendix 
C, we saw very complicate structure of more general results. Oscillation formulas typically involve all possible low- 
frequency and high-frequency combination terms. The amplitudes of the oscillation terms are essentially momentum 
dependent. We have also discussed the existence of the coherent antiparticle beam generated from the starting 
definite-flavor particle beam and presented its dynamics. 

Our general approach does not require to use any specific continuous parametrization of the mixing group but 
directly takes the values of matrix elements. This allows an analysis to be carried out in a unified closed form as 
shown in Sections II and III. In general, it may be preferable to solve the mixing problems without going through the 
intermediate parametrization step for the mixing matrix. Even if one wants to use a specific parametrization scheme 
for the mixing matrix, it is rather straightforward to formulate our general framework into a symbolic calculation 
system, like maple or mathematica, and carry out extensive calculations involving mixing parameters in short period 
of time. Examples of such calculations are shown in Appendix C. 

The physical application of the above formalism can be seen in investigating the neutrino mixing, mixing of gauge 
vector bosons governed by the Weinberg angle in the electroweak theory as well as vector mesons such as the p and ui. 
It seems also possible to apply these results to consider nonperturbative quark-mixing effects in the Standard Model 
and provide partial summation of the regular perturbation theory in mixing degrees of freedom. For this purpose 
considering covariant form of the above theory might be of great interest. Consideration along this line is in progress. 



14 



ACKNOWLEDGMENTS 



We thank Prof. G.Vitiello and Dr. M.Blasone for helpful discussions. This work was supported by a grant from the 
U.S. Department of Energy (DE-FG02-96ER 40947). The North Carolina Supercomputing Center and the National 
Energy Research Scientific Computer Center are also acknowledged for the grant of computing time. 

APPENDIX A: ESSENTIAL CASES OF MIXING FIELD PARAMETERS 

The most essential cases in modern particle physics are scalar/pseudoscalar (spin 0), vector (spin 1) boson fields 
and spin h fermion fields. For these cases mixing theory parameters are explicitly derived from quantum field theory 
[§,|l|. We then have for scalar/pseudoscalar fields (spin 0): 

u k,a = v k,o = !> ( A -!) 

for vector fields (spin 1): 

( fe ■elk) -\ 

(A.2) 



J k,0 \m> 1 rn "'J ' 
U k,±l = V k,±l = (°^«±): 



where n = ^ = e 



and ft± = ^f-^=(e x ± ie y ) form a spherical basis. For bispinor fields (spin 1/2), we use the standard 
representation of the 7-matrices given by 

(A.3) 



(A.4) 



. -/ J " V 
and the corresponding representations of spinors 



v -k.<j = ( _ \A( fc ) -m(n- g)uj- a , y / e(k) + muj- a ), 
where u>„ is spinor satisfying {n ■ a) Lo a — a ■ u a and a takes values ±1. 
The H and h matrix parameters are then for scalar case: 



for spin 1: 



h m _ h M = 1; (A.5) 



rrflj _ e A ,(fc)ej(fc)-fc 

k,0 m^vrij ' _ „ 

UM _ e M (fc)ej(fc)+fc 2 ' q\ 

k.O m^mj 



and for spin 1/2: 



= x /(e M (fc)+m At )(e J -(fc) + mj) + ^(e^ (fc) - m M ) (e,- (£:) - mj ), 



^ = o- (VOp W-™).)feW+ ffl j) - v 7 ^ (*0 +«v) (ej (fc) - m/jj 



(A.7) 
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APPENDIX B: THE FLAVOR VACUUM STATE 



In this appendix we explicitly solve flavor vacuum structure. We first consider boson case. 

We write the sought flavor vacuum state as the most general linear combination from the original-field Fock space 

1 



|0> 



r- E ■ 

(»),(*) 



— B(nm ( t) \..(4) h (buy...(bi k y\o 



(B.l) 



From the particle/ antiparticle symmetry, the part of Eq.(2.31) involving antiparticle annihilation operators results in 



a dependent set of equations and thus can be omitted. Expanding Eq.(2.31), we find: 

E K B (n 3 +l)0) + Pij B (n)(lj-l)) = °> a11 i n ) j (0 



(B.2) 



where (rij + 1) notation stands for (ni,n2, . . . ,rij + 1, . . .nfc) and fc is the number of flavor fields. To solve this 
infinite set of equations we introduce symbolic operators which decrease the subscript index of B coefficients, i.e. 



d-jBr n \(j-j — Br n \(j._iy Then solving each set of equation in (B.2) with respect to B ( n +i)(i} we nn( i 



-B (ni+ i)(2) = {J2 z ij d -j) B (n)(i) and consequently 



B (n)(l) "IKE Z H d - 



<B 



(0)0) 



(B.3) 



with matrix Z — — dr 1 ■ (3. Considering the momentum conservation and the original Eq.(B.2), it can be shown that 
only S(o)r;=o) must be non-zero among all (I). Thus, applying symbolic operators d-j and leaving only terms -B(o)(o) 
in the expansion, we get 

nA 



B 



(n)(Z) 




Z tk B (0)(0)- 



(B.4) 



It is possible to rewrite this complicate expression in the more compact form; 

1° >f= \ E II p(E Z »4blj) ki |0 >m (B.5) 
(fc) i j 

that can be shown directly by expanding the above expression. It also can be argued that to obtain Br n \m from 
Eq.( p3.^ ) one needs to leave only those terms in the expansion that give correct power of particle and antiparticle 
creation operators, i.e. total powers of all aj's are rij's and &t's are But this is same with extracting S(„vn 
from Eq.( |B.3] ). The constant Z is introduced instead of B(o)(o) and serves as a normalization factor determined by 
/<0|0>/=1. 



The Eq.(BJ5) can be further simplified as 



10 >/ = i EH ^T(E Zyalbi^lO > m = 
(fe) i j 

= in E ^r(E^«l^-) fc< |o> m = 

i k i= j 
N 

= ^exp( E Zijolbl^lO > m . 



(B.6) 



1G 



Let us now proceed to the fermion case. We employ the same idea with the symbolic shifting operators. If (?(„)(;) 
stands for creation operator for fermion state \(n), (I) >, we want then 



a i-S(n, + l)(Z)C(n, + l)(0 |0 >m— ±B( n . +1 ) (;) C(„) (j) 1 > m = d|_,-B( n ) (;) C(„) (j) |0 > r 
^^(n)(ii-l)A™)(^-l)l^ > ™ = ±- B (n)(i,-l)C'(n)(0|0 >m= d-iB( n ) (;) C*( n ) (;) 1 > T 



(B.7) 



with correct sign. Eq.(B.2) then can be written in the form 

J2( a H d +j + M-j) B {n)[l) = (B.8) 

3 

which binds together the shifting operators that increase and decrease the index. This set can be solved as 

d+i[B( n )(l)] = Z iJ d -A B (n)(l)] ( B - 9 ) 

3 

with the same matrix Z presented in the boson case. From the definition of shifting operators it can be inferred 
that they obey anticommutation property (i.e. d±id±j — —d±jd±i) and thus it can be shown further that for 
il > i 2 > ■ ■ ■ > i n 

d + i n d +in _ 1 . . . d+i 1 S( )( i ) = / B 1q n 

d-i t d-i 2 . . .«f_j,.B( n )(x) = B(„)(;_i) 

so that the solution can be written again as 

i 3 

where only i3( )(o) survives. Here, rii can be only or 1 and the anticommutation rules for the ordering are applied. 



It is remarkable that Eq.( B.5 ) can still be used for the fermion vacuum. This can be verified by a direct expansion 
with the anticommutation nature of ladder operators. Thus, for either boson or fermion case the flavor vacuum state 
can be written as 

1 N 

|0 >/= - exp( J2 ZiAb-jW >m ■ (B.12) 
We now proceed to find the normalization constant Z. To do this, we consider 

N N 

||0 > f | 2 = \eMj2 ZiAbljW >m I 2 = E IjaKE ^4^)^10 > m | 2 , (B.13) 

i,j = l L ' i,j=l 

where we use the fact that the states of (X^f^i ^ij a l^-j) L \^ >m are orthogonal for different L's. We then employ 
the fact that matrix Z can be transformed to a diagonal form with two unitary transformations, i.e. 

(B.14) 
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N 

We can now introduce additional unitary transformations of a'* = WaJ, b'> = to make Zija\b_- = 

HJ=1 

N 

^2 Z' ii a i ; where 6' satisfy the standard commutation/anticommutation relationship. Then, using the binomial 

i=i 

N 4. J. 

7/ 'ti,'t \L 



formula to expand Z^a^b_A , we find 



i=l 

L i=l 

-ExF E i! 2 n^l(^jO?6^) B, l >«l a = ( B -!5) 

L ni + ...+njv— L j — 1 

AT , 2 

= e e n^\z';/\ 2 = e 

L 'ni + ...+njv— L j — 1 ni,...,n^v 



where Aj's are eigenvalues of ZZ* . The summation limits in Eq.(B.lJi) are different for fermions and bosons. For 
bosons rii runs from to oo, while for fermions they only can be or 1. In either case the sum can be evaluated to 
give 

N ( IKH-Ai) fermions , ~ 77tw 

APPENDIX C: UNITARY MIXING OF 3 FIELDS IN WOLFENSTEIN PARAMETRIZATION 

We now present application of the above general formalism to specific case of mixing of 3 quantum fields. Calcula- 
tions were carried out with the help of Mathematica 3 symbolic calculational system. 

We note that all time-dependent quantities in this section are presented in the form of matrices each entry of which 
corresponds to certain flij = u>i + ujj or u>ij — LUi — LOj frequency. It means that each quantity is presented in the form 

(C.l) 




where P n and P w matrices are written as follows: 

nil _ r r pf! pf! t p£2 pf! pfii r pfi pfi p!! 11 

^ I \ 1 1 ■ 12 ■ 13 J ' I 21 5 22 ' 23 J ' I 31 ' 32 ' 33 J J 

P" = {{ift.ift, {^21^22, {ift, (C2) 

Since the diagonal elements of P w corresponds to the same constant term ua = 0, we can collect the diagonal elements 
of P u as Sp (P u ) = Pii + P22 + P33 and express only the off-diagonal elements as 

P" = {{0, PS, Pr 3 }, {P£, 0, P 2 - 3 }, {P3I, P31, 0}}. (C.3) 
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1. The case of 3 fermion fields 



We now show our results for the unitary mixing of 3 fields with spin i(bispinors). Although an explicit parametriza- 

tion is not needed in our formalism, we may use Wolfenstein parametrization as an explicit form of a mixing matrix 

/ 1 - A 2 /2 A AX 3 (p - in) \ 

U = -A 1 - A 2 /2 AX 2 . (C.4) 

\ AX 3 (1 - p - in) -AX 2 1 J 

All results are then computed to a few lowest orders in A. 

For the bispinors, we redefine our H and h matrices as — > / (2^/eTej), h 1 ^ — > (2^/eTej) so that 

(1 U12 Ul3 
"12 1 "23 

"13 "23 1 

Ui2 "13 
-"12 W 2 3 

~"13 -"23 



(C.5) 



Also, iXi j j "^ij tire defined in the same way as in the 2 field mixing 

■ v /(et+m i )(e J -+mj)+^/(e i -m i )(ej-mj) 



_ ^■ v /(e i -m i )(£j+mj)— v /(e i +mi)(ej-mj) 



(C.6) 



Then, the structure of the ladder operators is described by a and (3 matrices 



1 - A 2 /2 "i 2 A "13AA 3 (p - in) 

a = ( -U12A 1 - A 2 /2 U23AX 2 

u 13 AX 3 (1 - p - in) -u 23 AX 2 1 

/ " 12 A v 13 AX 3 (p-in)\ [ ' 

13= "12A " 23 A\ 2 

\ -"13AA 3 (1 - p - in) V23AX 2 

To make the results more compact, we define c = A(p — in), e = —A (1 — p — in) and a = A so that 

1 - A 2 /2 u 12 X "i 3 cA 3 

a = ( -"12A 1 - A 2 /2 U23C1A 2 

-ui3eA 3 — u 23 aX 2 1 

W12A "13CA 3 
(3 = \ v 12 X "23aA 2 

"i3eA 3 V23C1X 2 



(C.8) 



For the case when #2 flavor particle was initially present, the flavor charge oscillation formulas are as follows. The 
flavor charge fluctuation, Q2i2(t), is given by 



Q212 - {{0, (1 - A2 ) - K3"12 + "12"2 3 )}, 

{-^ (1 - A 2 ) , 0, ("12"13C* - "13"12C)}, (C9) 

( U23Vu + Ul2U23 ) i ( Ul2Wl3C * _ u 13 V 12 c) , 0}}, 



Q212 - {{0, (1 - A 2 ) , (-"12"23 + "12"2 3 )}, 

A 2 tt? 



j AfoMiacl (_ Ml2M23 + „ 12U23 ) ; ( Ml2Ul3C * + Vl2Vl3 c) , 0}}, 

^(Q 2 J i2) = A 2 (^2+^ 2 )(l-A 2 ). 
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Similarly, Q 222 (t) and Q 2 32(t) are given by 

Qg 2 = {{o, (i _ A 2 ) , " 2 (^ 3 ^2+ Ml2 ^3) 2 A 6 }i 

{ ^(l_A 2 ),0,^|^(l-A 2 )}, (C.11) 

r a 2 (M23"12+"12f23) 2 ^6 a 2 ^23 A4 ^ — A 2 ) 0}} 



Q% 22 = {{0, Hlli. (1 - A 2 ) , ° ("'3""-"""'3) A 6 }; 

{^(1- A 2 ), 0,^(1 -A 2 )}, (CJ2) 

f a 2 (tt23"12-f 12f23) 2 jy6 a "23 A ^ — A 2 ) 0}} 

Sp (Q& 2 ) = I - A 2 + i (3 '+ 2uf 2 + 4u 2 2 « 2 2 '+ 2« 4 2 ) A 4 , 



and 



Q232 = {{0, (e«i 3 «23 - w 23 « 13 e*) , ( U23 „ 12 + Ul2 „ 23 )} 7 

(eu 13 v 23 - u 23 v 13 e*) , 0, -^f^ + ^1}, (C.13) 

{ (U23U12 + Ui2W23) j + ii!^, o}}, 



Q232 = {{0, (ev 13 v 23 + u 13 u 23 e*) , (u 12 u 23 - «i 2 « 23 )}, 



( e *v 13 v 23 + u 13 u 23 e) , 0, -^f±- + £L^ }j 

,6 2 2 ,4 2 2 ,6 l^- 1 V 

{ Af^ue (ui2U23 _ Wi2U23) ^ _ ou^A. + iLi^, 0}}, 

(Q&2) = « 2 (t& + *4) A 4 (1 - A 2 /2) , 

respectively. 

In more details, the dynamics is given by the following quantities. The non-equal time anticommutators are given 
by: 

Fu(t) = e~ leit + A 2 (- e -* £lt + u 2 12 e~ le2t + v\ 2 e lf -^) , 
F 12 (t) = F 21 (t) = Xu 12 (e- l£2t - e-" 1 *) + A 3 i|a (e- 4£lt - e^ 2 *) , 
Fi 3 (t) = F 31 {-t)* = A 3 (w 13 (ce-" 3 * - e*e- tei *) - aui 2 u 23 e- 4e2 ' + av 12 v 23 e te2t -) , , 

F 22 {t) = e-" 2 * + A 2 (-e- l£2t + u 2 ^-^ + v 2 2 e* eit ) , 1 1 

F 23 (t) = F 32 (t) = X 2 au 23 (e-^ - e"^ 4 ) , 
F 33 = e— 3 *; 



Gii(t)=A 12l ; 12 (e-^-e^*), 
G12W - -(G 2 i(t))* = Avi 2 (e" iei ' - e" 2 ') + A 3 ^ 2 (-e^ 4 + e" 2 ') , 
Gi 3 (t) - -(Gsi(*))* = A 3 («i 3 (ee— 14 - c*e ie »*) + auss^e* 3 * + aw^^se"^ 2 *) , , p 1 

= A 2 Ml2Wl2 (-e"*«« + e" 1 ') , 
G 23 (*) - -(G32(t))* - A 2 ai; 2 3 (e~" 2t - e" 3 *) , 
G 33 (t) = A 4 a 2 U2 3«23 (e K2t - e^ 2 ') . 

The vacuum structure is defined by Z matrix: 

Zn = u 12 v 12 X 2 + U12 (l - u\ 2 ) wi 2 A 4 , 
Z12 = -v 12 \ - (| - m 2 2 ) wi 2 A 3 = Z 2 i, 
Z 13 = - (cv 13 - au 12 v 23 ) A 3 , Z31 = - {au 23 v 12 + ev 13 ) A 3 , 

Z 22 = -Zn + a 2 u 23 v 23 \ 4 , (C17) 
^23 = -aw 23 A 2 - (cu 12 v 13 + a (5 - u 2 2 ) w 23 ) A 4 , 
^32 = -aw 23 A 2 - (etii 3 + au i2 u 23 ) «i 2 A 4 , 
Z 33 = -a 2 U23«23A 4 . 

The normalization constant is obtained as Z rts 1 + v 2 2 A 2 + (y\ 2 + a 2 v 23 — v 2 2 u 2 2 )\ 4 + . . .. 
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If the particle of sort #2 was originally present, then for the particle of sort #1 the mixing quantities are as follows. 
The free-field particle condensate is 



Z[ = v 2 12 X 2 (C.18) 

and the flavor particle condensate is 

Z? = {{0, -A 2 4(l - A 2 ), -^v 13 c* (au 23 v 12 + v 13 e*)}, 
{-A 2 4(l - A 2 ), -X 4 uj 2 v 2 12 , -^aa (cu 13 v 12 + u 12 v 13 c*)}, (CJ9) 
{-^2-v 13 c* (au 23 v 12 + v 13 e*) , ~^-av 23 (cu 13 v 12 + u 12 v 13 c*) , 
-a 2 cu 13 u 23 v 13 v 23 c* A 10 }} , 



2^ = ||Q neui2^i3^23 \6 acMi 3 iii2«23 \6 j r ae"u 12 v 13^23 \6 Q 

a8 " 12 " 13C ' (2CU12U13 - 5^23) } ; { ac *" 13 2 " 12 '" 23 A 6 , 



2 A^psc (2c*u l2 u l3 - \au 23 ) ,0}}; 
^(Zn-A 2 « 2 2 -(1- U 2 2 )« 2 2 A 4 . 

The flavor particle number fluctuations are given by JV 2 i 2 (t) = \F 12 (t)\ z + Z^t): 



X a 



N? 12 = {{-\ 4 u 2 2 v 2 12 ,-\ 2V -?(l-\ 2 ), 



N 2 i 2 — {{0, -\ 2 u\ 2 (l — A 2 ) , 
V 1 (wi2«23 (—"i 3 e* + u 13 c - 2au 12 u 23 ) - cu 12 u 13 w 23 )}, {-A 2 w? 2 (l - A 2 ) 
0, {cu 12 u 13 + c*v 12 v 13 )}, 



(C.20) 



2 O13V23W12 (e* - c*) - v 12 v 13 u 23 c* + 2av 12 v 23 m 2 u 23 )}, 

{-A 2 4(1-A 2 ),0,-^A 4 }, (C.21) 
(*>13«23W12 (e* - c*) - W12W13W23C* + 2awi 2 w 23 ui 2 u 23 ) , 
-^AVaS^A 8 }}; 



(C.22) 



(wi2«23 (-u 13 e + u 13 c* - 2au 12 u 23 ) - c*u 12 u 13 u 23 ) , 
A!^ ( c *u 12 u 13 + cv 12 v 13 ),0}} ; 
Sp (JVft 2 ) = A 2 (« 2 2 + « 2 2 ) + A 4 (a 2 ^ - v\ 2 - u\ 2 (1 - v 2 2 )) , 

and the flavor antiparticle number fluctuations, N 2 i 2 (t), are given by 

AjQ _ f r ,,2 ,.2 \4 A 6 aeMi3t>i 2 f23 A 6 aMi 2 t> 2 3 /O^,, „ 1 a *\\ 
-< V 212 — II U 12 W 12 A ) 2 ' 2 \ Zau 23 v 12 + v 13 e )}, 

f A 6 aeu 13 v 12 V23 Q g t, 23 A 1 (C.23) 

{ A^a^ (2a U23 « 12 + v 13 e*) , -^1, -aSf^A 8 }}; 
^212 = {{0, - ^v^v 23 x^ _ av 12 ^ {2aUl2 u 23 + u 13 e*)} 7 

{- ae, " 12 ™ A6 ,0, 0}, {- a ™ A6 (2a Ul2U23 + ui 3 e) , 0, 0}}; (C24) 
(JV? 12 ) = A 4 (u 2 2 v 2 12 + a 2 v 2 23 ) . 

In the same initial condition, wc obtain the following for the particle of sort #2. The free-field particle condensate 

is 

Z' 2 = v 2 2 \ 2 +a 2 v 2 23 \ i (C.25) 

and the flavor particle condensate is 
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^ = {{-^ 2 2 A 4 ,-^(l-A 2 ) 
±f {2au 12 u 23 v 12 v 23 - c*u 23 v 12 v 13 + e* u 12 v 13 v 23 )} , 

„2„,2 



{^f (2au 12 u 23 v 12 v 23 - c*u 23 v 12 v 13 + e*u 12 v 13 v 23 ) , 
^X\-a^ 23 v^}}; 



{_^(1_A 2 ),0,-^3A 4 }, (C.26) 



Z$ = {{0, _ aeu 12 v^v M ^ j o^yAf _ 2a u 12 U 23 - e*U 13 )}, 

f ae*Mi2l)13f23^ 6 Q aU23V 12V 13C*A 6 j 

{^f^ (c*ui 3 - 2au 12 u 23 - eui 3 ) , ™ia2^cA! ; }}; 
5p (Z%) = v\ 2 \ 2 + {a 2 v 2 23 + v\ 2 (u\ 2 - 1)) A 4 . 



(C.27) 



The flavor particle number fluctuations, N 222 (t) = \F 22 (t)\ 2 + Z 2 (t), are given by 

= {{-eu 12 u 13 v 12 v 13 e*\ s , g™ A 6 , 
(° ("23^12 + u 12 v 23 f - u 23 v 12 v 13 c* + u 12 v 13 v 23 e*"j }, 

| gi3p3 A 6 |0| _^3 (cm 3 «i2 + (C28) 
(a (u 23 v 12 + u 12 v 23 f - u 23 v 12 v 13 c* + ui 2 vi 3 v 23 e*^j , 
_A^23 ( CWl3 „ 12 + saa) ,-a 2 cui3«23Wi3«23C*A 10 }}; 

% 22 -{{0,^(l-A 2 ), 

^ (a (u 12 U 23 " «12«23) 2 + W13W12W23C - Wl3«12«23e*) }, 

{^(1- A 2 ), 0,^f^}, 

{^T \ a ("12^23 - W12V23) 2 + U 13«12W23C - u 13 v 12 v 23 e*\ , (C.29) 

99.4 

-,o}} 



Sp{N^ 2 ) = \+\vl 2 -l) A 2 + 

(f + 4 1 - < + ^ 2 2 + 4 + «s 2 3 ) a 4 , 

and the flavor antiparticle number fluctuations, N 222 (t), are given by 

^222 = {{-eui 2 ui 3 v 12 v 13 e*X 8 ,-^- (l - A 2 ) , ^v 13 (-w 23 vi 2 c* + wi 2 v 23 e*)}, 

(1 _ A 2) , 0j _^! } , { ^ Vl3 (~u 23 v 12 c* + u 12 v 23 e*) , (C30) 
-a 2 cc*ui 3 ii 2 3«i3W23A 10 }}; 

^222 = {{0, ~ ae " 12 ^ 13 " 23 A 6 , Q "13^12^3A 6 ( c _ e *)} ; 
/ ae*Mi 2 fl3f23 \6 n a "23 - »12'»13c'' A 6 1 

i 2; A > u > 2_ 7 1 ' (C31) 

|Otti3_Ui2t)23A_ _ ^ atl23"12" 13cA 0}}' 

Again in the same initial condition, the mixing quantities for the particle of sort #3 are as follows. The free-field 
particle condensate is 

Z[ = a 2 v 2 3 X 4 (C.32) 

and the flavor particle condensate is 



-^232 = {{-eui2Ui 3 vi 2 vi 3 e* A 8 , — (u 13 v 23 e + U23«i3e*) 
(au l2 v 23 - v 13 c*)}, {^12 ( Ul3 v 23 e + u 23 v 13 e*) , 



a 2 Ml2W 2 3«12W23A 6 , ^22.} 



(C.33) 



(au 12 v 23 - Wl3 c*) , -^J^,0}} 
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%32 = {{0,a 2 eu 13 u 23 v 13 V2 3 e*X w ,-^^f^X 6 }, 
' {a 2 eu 13 u 23 v 13 v 23 e*X 10 , , -^^i3£l A 6 }, 



2 

f au 13 v 1 2V23e \6 au 2 3V 12^13C \6 Q||- 

5p(Z 2 - 32 ') = a 2 4A 4 . 



(C.34) 



The flavor particle number fluctuations, N 232 (t) = \F 32 (t)\ 2 + Z 3 (t), are given by 

Ng 2 = {{-\*ul 2 vl 2 ,-\^(l-\ 2 ), 
±f (2av 12 v 23 u 12 u 23 - u 23 v 12 v 13 c* + v 13 {u 23 v 12 + u 12 v 23 ) e*)}, 

{-A 2 4(1-A 2 ),0,-^A 4 }, (C.35) 
(2avi 2 v 23 ui 2 u 23 - u 23 v 12 vi 3 c* + vi 3 {u 23 v 12 + u 12 v 23 ) e*) , 
^^x\-aW 23 v* 23 \Z}}; 

%S2 - {{0, (v 13 v 23 e + u 13 u 23 e*) , 
( w i2 w 23 ("13 (c - e*) - 2au 12 u 23 ) + u 12 u 13 u 23 e*)}, 
i-^ 1 (u 13 u 23 e + v 13 v 23 e*) , 

0,-^f^}, (C36) 
{^r («12«23 (wi3 (c* - e) - 2au 12 u 23 ) + ui 2 u 13 u 23 e) , 

2 2 \ 4 

_a^3A_ )0}} . 

(N% 32 ) = A 2 « 2 2 + A 4 (a 2 {vl 3 + u 2 3 ) - v\ 2 (l - M 2 2 )) . 
Similarly, the flavor antiparticle number fluctuations, N 232 (t), are given by 

JV?32 = {{-^2«? 2 A 4 , (1 - A 2 ) , A!£I^3 ( 2aW23Ul2 - v 13 c*)}, 

{-^f. (1-A 2 ) (C.37) 
{ A6a T M23 (2a«23«i2 - v 13 c*) ,- °™2^ ,-a 4 n 2 3 ^ 3 A 8 }}; 

^232 = {{0, 0, ( 2aUl2U23 - u 13 c)}, 

/A n aM 23« 12« 13C* A 6 1 

{0,0, 2 j, (C38) 

{- (2aW 12U23 - U 13 C*) , a M23 ^l3cA 

(7V 2 - 32 ) = « 2 2 A 2 + AS 2 2 (uj 2 - 1) . 



2. The case 3 boson fields 



We now consider the application to bosons. The boson case is not much different from the fermion case. With the 
use of the 7^,7^ matrices, one can write the ladder mixing matrices as 



l-A 2 /2 7 + A j+AX^p-ir,) 

a = ( -7+A 1 - A 2 /2 7 2 + 3 ^A 2 

7l + 3 AA 3 (1 - p - ir,) -7 2 V1A 2 1 

7i 2 A i 13 A\ 3 (p - irj) 

(3=\ 7l2 A 7 23^A 2 

--/ 13 AX 3 (1 - p - t V ) j 23 AX 2 



(C.39) 



We see that a and (3 indeed have the same form as in the fermion case with the correspondence 7^ — > Uij and 



7»j 
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1 - A 2 /2 u i2 A u 13 c\ 3 
a = I — U12A 1 - A 2 /2 ii 2 3oA 2 

"■■:< A 3 »-'.;"A 2 1 , r , 
' v 12 X v 13 cX^ 
/3 = v 12 X 

y «i3eA 3 w 2 3«A 2 

This shows that the only difference appears in the quantities that have explicit spin dependence, i.e. F^ v and 
everything involving F^ v . As a rule, the quantities for the boson case can be obtained from the fermion formulae by 
simply changing the signs in the terms quadratic in Vij . We summarize them below. 
The oscillation formulas are as follows. For the particle of sort #1, Q 2 i 2 (t) is 

Q212 = {{0, ^ (1 - A 2 ) , (u 23 V 12 + u 12 v 23 )}, 

(1 _ A 2 ) , 0, (u 12 v 13 c* - u 13 v 12 c)}, (C41) 

{ A6a 2 13C * («23«i2 + u 12 v 23 ) , -^23 ( Ul2 „ 13C * _ Ul3Wl2 c) , 0}}; 



Q212 = {{0, (1 - A 2 ) , -Aispc (ui2U23 + Ul2U23 ) }j 

{ _ ^ _ A 2) ^ Qj ( Ul2Ul3C _ Ul2Ul3C *)}, 

{- A "™" c ' ( M 12«23 + V 12 V 23 ) , (Ml 2 M 13 C* - V 12 V 13 c) , 0}}; 

^(Q 2 J 12) = A 2 («?2-^2)(1-A 2 ). 



r a (U23V12+U 12 V23) ^6 g "23 A fa _ ^2^ Qjj. 



0222 = {{0, (1 - A 2 ) , ° ("»3«12+«12«>23) A 6} ; 

(1 _ a 2 ) ,0, (! - A 2 )}, 

(- a 2 (tt23fl2+fl2f23) 2 ^6 ° 2 "23 A4 ^ — A 2 ) 0}}" 

Sp (Q% 22 ) = \ - A 2 + \ (3 '+ 2uf 2 - 4u 2 2 « 2 2 V 2Jf 2 ) A 4 . 



For the particle of sort #3, Q 232 (fy is 



The non-equal time commutators are given by: 



F n (t) = e- teit + A 2 (-e- leit + u\ 2 e~ lf ^ 



?if-2t\ 



(C.42) 



For the particle of sort #2, Q 222 (t) is 

Q222 = {{0, (1 - A 2 ) , - a2( " 23 " 12 2 +ttl2 " 23)2 A 6 }, 

{-^(l- A 2 ), 0,-^4^(1 -A 2 )}, (C43) 



(C.44) 



Q232 = {{0, (eu 13 v 23 - u 23 v 13 e*) ( U23 „ 12 + u 12 v 23 )}, 

{-^ (e Ul 3«23 - u 23 v 13 e*) , 0, - ^f^}, (C.45) 

Q232 = {{0, (ev 13 v 23 - u 13 u 23 e*) , A!s|3el + Vl2V23 )} i 

{ au^L ( e *v 13 v 23 - u 13 u 23 e) , 0, + g!^}, 

{ (ui2U23 + Wi2U23) ; + £i!43^, }}; 

Sp (Q& 2 ) = a 2 (u 2 23 - v 2 23 ) A 4 (1 - A 2 /2) . 



(C.46) 



F 12 {t) = F 21 {t) = \u 12 (e~ le2t - e- leit ) + A 3 ^f (e^ 1 - e^ 2 ') 
Fi 3 (t) = (F31 (-*))* = A 3 (cuise-^* - aUl2 u 2 3e- je2t - av^s^ 2 * - e* Ul3 e-* eit ) , , , 

F22 (t) = e^ £2t + A 2 (-e-* £2t + u 2 2 e- jeit - v 2 2 e lf -^) , 1 j 

F 23 (i) = F 32 (t) = X 2 au 23 (e-"3* _ e -* 2 t) , 
*33 (t) - e-"3* ; 
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G n (t) = \ 2 u 12 v 12 ( e - fe * - e ie2t ) , 
Gv»{t) = -{G 2 i{t))* = Xv 12 (e- leit - e" 2 *) + A 3 ^ (-e" 161 ' + e l£2 *) , 
Gi 3 (t) = -(G 3 i(t))* = A 3 (cm 23 t>i 2 e I£2t + ewig • e" i£li + auia^e"" 2 * - u 13 c*e^*) , . . 

G 22 (t) = X 2 u 12 v 12 (-e-^ + e^*) , l ° 8j 

G 23 (t) = -(G 3a (t))* = ^23 (e" l£2t - e <e »*) , 
G 33 (t) = A 4 a 2 U23U23 (e^'-e— 2 *). 

The vacuum structure is given by the fermion Z (Eq.( |C.i"7| )) with the normalization constant 

-3-^11 2 \ 2 , 2 . 22 i 4 2 2 \ \ 4 , 

Z « 1 + v 12 X + {v 12 + a v 23 + v 12 - v 12 u 12 )X + ... 

If the particle of sort #2 was emitted initially, then the particle of sort #1 has the following free-field particle 
condensate: 

Z[ = v 2 2 X 2 (C.49) 



and the flavor particle condensate identical to the fermion case, i.e. Eqs.(C.19) and (C.2C) 
The flavor particle number fluctuation N 2 \ 2 (t) is given by 



Ng 2 = {{-X^ 12 -\^-f{l-\^ 



^ (2au 12 u 23 v 12 v 23 + ui 2 v 13 v 23 (e* + c*) - u 23 wi2«i 3 c*)}, 



{_A 2 ^(l-A 2 ),0,-^§^}, (C.50) 



{^f (2au 12 u 23 v 12 v 23 + u 12 v 13 v 23 (e* + c*) - u 23 v 12 v 13 c*) , 

2 2 \4 



2 

-a 4 u 2 23 v 2 3 X 8 }}; 



% 12 = {{0,-A 2 u 2 2 (l-A 2 ), 
{v\ 2 v 23 (u 13 e* - u 13 c + 2au 12 u 23 ) + cu 12 u 13 u 23 )} , {-X 2 u 2 2 (l - A 2 ) 
0, (cu 12 u 13 + c*v 12 v 13 )}, 

{-^f - (vi 2 v 23 {u 13 e - u 13 c* + 2au 12 u 23 ) + c*u 12 u 13 u 23 ) , 
(c*u 12 u 13 + cv 12 v 13 ) , 0}}; 
Sp (JVft a ) = A 2 (u\ 2 + v 2 2 ) + A 4 (a 2 v 2 23 v\ 2 - u\ 2 (l - v\ 2 )) . 

The flavor antiparticle number fluctuation N 2 \ 2 {t) is given by 



{~2u 23 v 12 v 13 c* + u 12 v 23 (2au 23 v 12 + Ui 3 e*)) , 



(C.51) 



^12 = {{-" 2 2« 2 2A 4 ,-« 2 2 A 2 (1-A 2 ) 
(~2u 23 v 12 v 13 c* + u 12 v 23 (2au 23 v 12 + v 13 e*))}, 

{-v 2 2 X 2 (l-X 2 ),0,-^f^l (C52) 



#212 = {{0, - ^pA! r aiSA! (-2c Ul3+ 2a Ul2M23 + M i 3 e*)}, 

{- ae '" 12 " 2 13 " 23A6 ,0, a«23«i2«i 3 c*A 6 }, (c 53) 

{-^f^{-2c*u l3+ 2au 12 u 23 + u 13 e),au 23 v 12 v 13 cX\Q}}- 
Sp (NZ 12 ) = 2v 2 2 X 2 + A 4 (u 2 2 t; 2 2 + a 2 v 2 23 - 2v 2 2 ) . 

For the same initial condition, the particle of sort #2 has the free-field condensate given by: 

Z' 2 = v 2 2 X 2 + a 2 v 2 3 X 4 (C.54) 
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and the flavor particle condensate identical to Eqs.(C26) and (C.27). 
The flavor particle number fluctuation N 222 {t) is 

^ 22 = {{-2^ 2 A 4 ,- W ? 2 A 2 (l-A 2 ), 

(a (u 23 V 12 - M12W23) 2 + U 23 V 12 V 13 C* - Ui 2 V 13 V 23 e*^j }, 

{-« 2 2 A 2 (l-A 2 ),0,-a 2 t4A 4 }, (C55) 
{-^f ( a ( M 23«i2 - M12W23) 2 + u 23 vi 2 v 13 c* - u 12 vi 3 v 23 e*^j , 
- a 2 « 2 3 A 4 ,-2a 4 u 2 3U 2 3 A 8 }}; 

% 22 - {{0, ^ (1 - A 2 ) , 

(cMi 3 Wi2f 2 3 + a (U12U23 - Wl2f23) 2 ~ ^13^12^236*) }, 

{^(1_A 2 ),0,^K}> 
{^T { c *ui 3 vi 2 v 23 + a (U12U23 - W12W23) 2 - ui 3 vi 2 v 23 e \ , (C.56) 

2 2 \ 4 

^f^,0}}; 
Sp(N$ 22 )=% + (vf 2 -l) A 2 + 

! + 4* - Vh + U\ 2 V\ 2 + ^ + « 2 «2 2 3 ) A 4 - 

The flavor antiparticle number fluctuation N 222 {t) is 

^22 = {{-2^ 2 2 A 4 ,-^(l-A 2 ), 

V s (-"23«i2Wi 3 c* + U12W23 (4au 23 wi2 + v 13 e*))}, 

{ _^!(1_ A 2) j0j _^1 }; (C.57) 



(-W23W12W13C* + Mi2^23 (4aM 23 vi2 + «i3e*)) , - ° v f , 
-2a 4 U 2 3 <; 2 3 A 8 }}; 

N 222 = {{0, - ae " 12 ^ 23 A 6 , 221^ (( c - e *) Ul3 - 4au 12 u 23 )}, 

/ ae*u 12 v 13^23 \6 n g "23'"12fl3C* A 6 1 

J 2 A ' U ' 2 h , 6 (C58) 

{ a " 12 r ((c* - e) u 13 - 4 aUl2M2 3) , aU23Vl2 2 13cA , 0}}; 
5p (#& 2 ) = « 2 2 A 2 + (a 2 W 2 3 - v\ 2 + 1u\ 2 v\ 2 ) A 4 . 

Finally, for the same initial condtion, the particle of sort #3 has the free-field particles condensate given by 

Z[ = a 2 w 2 3 A 4 (C.59) 



and the flavor particle condensate identical to Eqs.( C33 ) and ( C.34 ). 
The flavor particle number fluctuation N 232 (t) is given by 

^32 = {{-AV 2 « 2 2,-A 2 4(l-A 2 ), 

^f- (~U 23 Vi 2 V 13 (c* + e*) + 2aUi 2 M 2 3Wl2«23 + e*Ui2fl3W 2 3)}, 

{-A 2 4 2 (l-A 2 ),0, 

(-M23V12U13 (c* + e*) + 2aui 2 w 23 Wi2W23 + e*Mi2Vi3«23) , 
-^A 4 ,-a 4 u 2 3t ; 2 3 A 8 }}; 

^232 = {{0, -22^1 (v 13 v 23 e + u 13 u 23 e*) , 

V s (^12^23 («13 (C - e*) - 2aUi 2 U 2 3) + Wl3Wl2«23e*)}, 

6 {-^f^ {v 13 v 23 e* + u 13 u 23 e) , 0, -^^+, 

{^T ( V 12V 23 (ttl3 (C* - e) - 2aUi 2 U23) + Wl3Wl2U23e) , 

2 2 \ 4 

-^f^,0}}; 

Sp (N& 2 ) = \ 2 v\ 2 + A 4 (a 2 (u 2 3 + W 2 3 ) - vl 2 (1 - u? 2 )) . 



(C.60) 



(C.61) 
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The flavor antiparticle number fluctuation N232 (t) is given by 



^ 2 = {H%^,-A 2 ^(1-A 2 ), 

(-U23W12W13C* + 2aUi 2 U23Wl2f23 + 2e*Ui 2 Wl 3 W23)}, 

{ _ A2 4(l_A 2 ),0,-a 2 U2 2 3 A 4 }, ( C - 62 ) 

(-U 23 V 12 VlsC* + 2aUi 2 U 2 3Wl2W23 + 2e*1ti 2 ^13 «23 ) , 

-a 2 w| 3 A 4 , -a A ul 3 vj 3 X 8 }}; 

^232 = {{0, -A 6 aui 2 wi 3 W23e, 
^^12^23 ("13 (c - 2e*) - 2aui 2 u 2 3)}, 

{-A 6 aUl2 «i3«23e*, 0, Alf^c*^^}, (C.63) 
{^12^23 («i3 (c* - 2e) - 2au 12 U23) , ^^cv 12 v 13) 0}}; 
^ (^32) = A 2 « 2 2 + A 4 (2a 2 V 2 3 - « 2 2 (1 - ul 2 )) . 
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